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Abstract 

We investigate the Affleck-Dine baryogenesis after D-term inflation with a positive 
Hubble-induced mass term for a B — L flat direction. It stays at a large field value 
during D-term inflation, and just after inflation ends it starts to oscillate around the 
origin of the potential due to the positive Hubble-induced mass term. The phase 
direction is kicked by higher-dimensional Kahler potentials to generate the B — L 
asymmetry. The scenario predicts nonzero baryonic isocurvature perturbations, which 
would be detected by future observations of CMB fluctuations. We also provide a 
D-term inflation model which naturally explain the coincidence of the energy density 
of baryon and dark matter. 


1 Introduction 


Inflation is a new paradigm to solve cosmological problems related to the initial conditions of 
the early Universe. However, any preexisting baryon asymmetry is diluted away due to the 
inflationary expansion of the Universe, so that there should be a mechanism to generate the 
observed baryon asymmetry after inflation. In addition, the observed abundance of baryon 
asymmetry is equal to that of dark matter (DM) within of order unity (H fe /D DM ~ 0.2) [1], 
which is a mystery in cosmology referred to as the baryon-DM coincidence problem. This 
implies that the baryon asymmetry and DM have a common origin. We therefore need to 
consider inflation, baryogenesis, and DM production simultaneously to construct consistent 
cosmological models. 

In this paper, we focus on D-term inflation models by the following reasons [2, 3]. First, 
the energy scale of D-term inflation is naturally of order the GUT scale, which predicts an 
amplitude of CMB fluctuations consistent with the observations [4, 5]. The second reason 
is related to the so-called 77 -problem. If inflation is driven by a nonzero F-term potential 
energy, supergravity effects induce masses of order the Hubble parameter to all scalar fields, 
including inflaton. However, the Hubble-induced mass for inflaton spoils the flatness of its 
potential and results in an 0(1) slow roll parameter 77 ~ 1 . Since Hubble-induced masses 
come only from nonzero F-term potential energy, the 77 problem is absence in the case of 
D-term inflation, in which inflation is driven by nonzero D-term potential energy. 

In supersymmetric (SUSY) theories, baryon asymmetry can be generated by the AfHeck- 
Dine baryogenesis using a B — L charged flat direction called an AD field [ 6 , 7]. Since the 
soft mass of the AD field is much smaller than the energy scale of inflation, it can stay at 
a large vacuum expectation value (VEV) during D-term inflation. After D-term inflation 
ends, the AD field obtains a Hubble-induced mass from the oscillating inflaton held through 
the supergravity effects. In the literature, the Hubble-induced mass term was assumed to be 
negative (i.e., tachyonic) or absent to investigate the dynamics of the AD held and calculate 
the produced baryon asymmetry [ 8 , 9, 10, 11]. In this paper, we investigate the Afheck- 
Dine baryogenesis in the case that the AD held obtains a positive Hubble-induced mass 
term after D-term inflation. In this case, the AD held starts to oscillate around the origin 
of the potential due to the positive Hubble-induced mass term just after inflation. At the 
same time, the phase direction of the AD held feels non-renormalizable B — L breaking 
operators and starts to rotate in the phase space, which results in a generation of B — L 
asymmetry. Then, the coherent oscillation of the AD held decays and dissipates into the 
thermal plasma and the B — L asymmetry is converted to the desired baryon asymmetry 
through the sphaleron effects [12, 13]. We calculate the baryon asymmetry and show that 
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the result can be consistent with the observed amount of baryon asymmetry. Since the 
AD field fluctuates due to the absence of Hubble-induced mass during D-terrn inflation, the 
Affleck-Dine baryogenesis predicts some amount of baryonic isocurvature perturbations [9, 
10, 11, 14, 15]. Especially, when the AD field starts to oscillate due to the positive Hubble- 
induced mass term, its radial direction also has quantum fluctuations and contributes to 
the baryonic isocurvature perurbations. If one might consider a high-scale D-term inflation 
model, the resulting isocurvature perturbations would be detected by future observations of 
CMB fluctuations. 

We also build a D-term inflation model which naturally predicts the baryon-to-DM ratio 
of order unity. We introduce a shift symmetry for the inflaton superfield to ensure a large 
initial VEV of inflaton [17]. We also introduce its linear term in the Kahler potential so that 
the branching of inflaton decay into gravitinos can be of order unity [18, 19, 20, 21, 22, 23, 24], 
When the mass of gravitino is larger than 0(100) TeV, it decays into the minimal SUSY 
standard model (MSSM) particles before the big bang nucleosynthesis (BBN) epoch. The 
DM, which is the lightest SUSY particle (LSP), is therefore produced non-thermally from the 
gravitino decay. This scenario, together with the above Affleck-Dine baryogenesis scenario, 
predicts an 0(1) ratio of the energy density of baryon and DM. This means that the scenario 
naturally explains the baryon-DM coincidence problem. This arises from the fact that both 
of them is related to the energy scale of inflation. The amount of baryon asymmetry is 
proportional to the reheating temperature of the Universe and inversely proportional to the 
Hubble parameter during inflation. The DM abundance is proportional to the reheating 
temperature and inversely proportional to the mass of inflaton. Since the Hubble parameter 
and the mass is related to each other, the resulting baryon and DM density is naturally of 
order unity. We predict that the LSP mass is two orders of magnitude larger than the proton 
mass, which comes from the fact that the GUT scale is two orders of magnitude less than 
the Planck scale. When the LSP is mostly wino or higgsino, it would be detected by future 
indirect detection experiments of DM. 

This paper is organized as follows. In the next section, we briefly review the simplest D- 
terrn inflation model as an illustration. In Sec. 3, we consider the Affleck-Dine baryogenesis 
and calculate the baryon asymmetry and baryonic isocurvature fluctuations in the case that 
the AD field obtains a positive Hubble-induced mass term after inflation. In Sec. 4, we 
provide a D-term inflation model which naturally predicts the baryon-to-DM ratio of order 
unity. Section 5 is devoted to the summary. 
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2 D-term inflation 


We focus on D-term inflation [2, 3], in which inflation is driven by a finite energy density of 
the D-term potential. Although the following simple model of D-term inflation predicts the 
spectral index relatively blue tilted compared with the observation of CMB fluctuations, we 
review it as an illustration. Note that there are variants of D-term inflation models which 
predict the spectral index consistent with the observed value [25, 26], and the results in the 
next section can be applied to those models, too. 1 

We introduce a £/( 1) gauge symmetry with a Fayet-Ilipoulos (FI) term £ and consider 
superfields S, and ip + with U( 1) gauge charges as 0, —1, and 1, respectively. 2 The 
D-term potential is written as 

Vb = y (foM 2 - foM 2 -0 2 > W 

where g is the 17(1) gauge coupling constant. We introduce a superpotential given as 

W (inf) = A Si/j+jJj-, (2) 


where A is a coupling constant. Hereafter, we denote their scalar components by the same 
symbols as the superfields. 

The scalar component of the field S plays a role of inflaton. Suppose that the inflaton S 
has a VEV larger than the critical value of S c = gy/£/\. The fields and rp + obtain large 
effective masses from the VEV of the inflaton and stays at the origin of the potential. In this 
regime, the nonzero D-term potential of Vo = g 2 £ 2 /2 drives inflation. The Coleman-Weinberg 
potential for the inflaton lifts its potential above the critical point such as 


Ei_ioop ^ ^g 2 i 2 



167r 2 


log 



( 3 ) 


where Q is a renormalization scale. Thus, the inflaton slowly rolls down to the origin of the 
potential. The COBE normalization requires [4, 5] 


>/% ~ 6.6 x 10 15 GeV. (4) 

1 The results are not applicable to the inflation model considered in Ref. [27] because a F-term potential 
drives inflation with a sizable e-folding number in that scenario. 

2 Thc FI term can be generated dynamically as discussed in Ref. [28]. They also argue that the production 
of cosmic strings at the end of inflation is generally avoided when the FI term is generated dynamically. 
Otherwise the CMB data puts an upper bound on a cosmic string contribution to the CMB fluctuations, 
which leads to the upper bound on the FI term [25, 29]. 
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This leads to the Hubble parameter during inflation such as 


H r ~ ^ ~ 3.7 x 10 12 GeV 


\/3M] 


pi 




6.6 x 10 15 GeV 


The e-folding number and a slow roll parameter are calculated as 


A* ~ 

V = 


4t r 2 S 2 
9 2 M&' 

V" „ t2 9 2 1 

V P1 ~ 8vr 2 SI 21V* ’ 


( 5 ) 

( 6 ) 
( 7 ) 


where the subscript * denotes values corresponding to the pivot scale &* = 0.05Mpc -1 . The 
slow roll condition fails (77 ~ 1) at the VEV around S ~ g/( 2 - 0270 , which is larger than the 
critical value S c for the case of A = 0(1). Thus, slow roll inflation ends at the VEV around 
S — g/(2\Z2ir) and soon after that the waterfall field 0 + starts to oscillate around the low 
energy minimum of i/£. 

The scalar spectral index is calculated as 

n s ~ ~ °- 98 - ( 8 ) 


It is measured by Planck data alone such as [5] 


n (° bs ) = 0.9616 ± 0.0094 ( 68 %). 


(9) 


So the above prediction deviates by about 2 cr from the observation. Let us emphasize that 
the results in the next section can be applied to other variants of D-terrn inflation models, 
including the ones which predict the spectral index consistent with the observed value within 
a la level [25, 26]. 

After inflation ends, the energy density of the Universe is dominated by that of the 
oscillation of S and ifj + . When some MSSM fields carry nonzero 17(1) charge, the field "0+ 
immediately decays into the MSSM fields through the interaction in the D-terrn potential. 
Even if the MSSM fields have no 17(1) charge, the kinetic mixings between the 17(1) and 
U(l)y makes the field 0+ decay into the MSSM fields relatively fast [ 8 ]. Thus, the reheating 
temperature of the Universe is determined by the relatively late-time decay of the inflaton S, 
which dilutes the relics produced from the decay of 0 + . We define the reheating temperature 
as 
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where Ts is the decay rate of the inflaton S. The reheating temperature of the Universe 
depends on the mass of the inflaton (ms = A\/£) and the assumption of interactions between 
S and the MSSM fields, which is determined by their underlying symmetry. We explicitly 
calculate the reheating temperature in Sec. 4 for a specific model with a shift symmetry and 
an approximate Z 2 symmetry, while we regard it as a free parameter in the next section to 
calculate the the amount of baryon asymmetry generated by the Affleck-Dine baryogenesis. 

3 Affleck-Dine baryogenesis 

In this section, we consider the Affleck-Dine baryogenesis [6, 7] after D-term inflation and 
calculate the resulting baryon asymmetry and baryonic isocurvature perturbations. We 
consider the case that the AD field obtains a positive Hubble-induced mass term after the 
end of inflation, which has been overlooked in the literature. We investigate the potential of 
the AD field after D-term inflation in the next subsection, and then we calculate the baryon 
asymmetry. In Sec. 3.3, we calculate baryonic isocurvature perturbations predicted by the 
Affleck-Dine baryogenesis. Finally, we comment on Q-ball formation in Sec. 3.4. 

3.1 Potential of the AD field 

We consider the Affleck-Dine baryogenesis [6, 7] using a flat direction (=AD field) (p with 
nonzero B — L charge. The AD field has soft SUSY breaking terms through the low-energy 
SUSY breaking effect. Since the soft mass of the flat direction is much smaller than the 
Hubble parameter during inflation, it has a large VEV during inflation. In this paper, we 
assume that the superpotential of the AD field is absent or sufficiently small so that the 
initial VEV of the AD field (pi can be as large as the Planck scale ((pi ~ M P i). 3 Such a large 
VEV is favoured to avoid the baryonic isocurvature constraint as shown in Sec. 3.3 (see also 
Ref. [15]). Note that owing to the exponential term in the supergravity potential the VEV of 
the AD field is restricted below the Planck scale. Since the curvature of the phase direction 
is absent (or at least much less than the Hubble parameter), the phase of the flat direction 
also stays at a certain phase during inflation. We denote the initial phase of the AD field as 

6i. 

3 For example, we can introduce R-symmetry in which the charge of the AD Held is zero to forbid the 
superpotential for the AD field. This symmetry is not exact because it is inconsistent with the constant term 
in the superpotential which is needed to ensure the (almost) vanishing cosmological constant as well as the 
gaugino mass terms. Since the R-symmetry breaking order parameter is of order the gravitino mass, which 
is much smaller than the energy scale of inflation, such breaking terms can be neglected in the following 
discussion. 


5 



After inflation, the energy density of the Universe is dominated by that of oscillating 
inflaton. Since the inflaton oscillation induces F-term potential, the flat direction obtains 
Hubble-induced terms through supergravity effects [7]. The scalar potential in supergravity 
is given as 


V = e K/M P' 


(ZW) K* (DjWy - J^ \W\ 


( 11 ) 


where K is a Kahler potential and DiW = Wi + K,W /Mp,. The subscripts represent the 
derivative with respect to field i and K u = (K,j)^ 1 . Since the F-term of is given by 
A Sijj +} the scalar potential includes the following term: 

I j 12 


V D 


Mk 


|A Sfa 


( 12 ) 


Using 0_|_ = y/£ and taking average with respect to time, we obtain the Hubble-induced mass 
term of 


V 3 |i/ 2 (l) \<t>\\ 


Here we have used the virial theorem: 


ml{\S\ 2 )~-H-(t)M 2 m , 


(13) 

(14) 


where () represents the time-average. In the literature, they assume a negative Hubble- 
induced mass term which comes from higher-dimensional Kahler potentials, say, 

K ‘ m = Cs |s|2 (15) 

where cs is an 0(1) constant. As we stated in the introduction, we consider the positive 
Hubble-induced mass term and do not need to introduce such higher-dimensional terms. 
Hereafter, we consider a general case and denote the coefficient of the Hubble-induced mass 
term as ch (= 0(1)). 

In addition to the Hubble-induced mass term, the flat direction obtains higher-dimensional 
terms from non-renormalizable Kahler potentials. 4 The following Kahler potential may exist 
and induce 17(1) breaking higher-dimensional potential for the flat direction: 






-a H / d 2 6d 2 9\S\ 


nMi 


+ c.c. 


pi 


—Off IcUSI 2 -1- c.c. 

3 1 M 1 nM, n 


pi 


-a H H 2 (t) 


nMpf 2 


+ c.c. 


( 16 ) 


4 The usual Hubble-induced A-terms are absent during the inflaton oscillation era. 
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where n is an integer depending on flat directions. For example, n — 3, 6, 9,... for u c d c d c 
flat direction. In the last line, we take average with respect to time and use the relation 
of ((<9o*S') 2 ) ~ 3H 2 (t)Mpi/2, which comes from the virial theorem. Note that this term has 
a nonzero phase which is different from the phase of the flat direction 6i during inflation. 
We can redefine the phase of the flat direction to eliminate the phase of a B - After the 
elimination, we redefine the initial phase of the flat direction as 9i without loss of generality. 
The discrepancy between the initial phase of the flat direction and the phase of the above 
1/(1) breaking term is essential to generate the baryon asymmetry. 

In summary, the AD field obtains the following potential after inflation: 

v(t) = c u H\t) M 2 - a„H'\t) + C.C.) + ..., (17) 

where c B and a B are positive 0(1) parameters. The dots represents higher-dimensional 
terms which restrict the AD field below the Planck scale. Since the flat direction starts to 
oscillate due to the Hubble-induced mass, we can neglect usual soft mass and A-terms for 
the AD field. 


3.2 Calculation of baryon asymmetry 


In this subsection, we calculate the baryon asymmetry generated from the AD field with the 
potential (17). The initial VEV and phase are (j) l and 9 i} respectively. For c H > 0, the flat 
direction starts to oscillate around the origin of the potential just after the end of inflation. 
At the same time, the flat direction is kicked in the phase direction due to the second term 
in Eq. (17). The B — L asymmetry is generated through this dynamics. The evolution of 
equation for the B — L charge density is written as 


Ab_l + 3 Hub-l 


2q Im 



(18) 


where q denotes the B — L charge of the AD field. From this equation we obtain 


a 3 n B -L(t osc ) ~ / dt2qa 3 (t) \(j>V' A \ sin (n9) 


= eqHjtf 
e ~ (3 — 4) x 


3n — 6 


a B sin (— n9i ) 


M, 


pi 


n —2 


(19) 

( 20 ) 


where we have used 0 oc a 3,/4 . We define e (< 1) which represents the efficiency of baryoge- 
nesis. We have numerically solved the equations of motion for 0 and S with the Friedmann 
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equation and have obtained the numerical factor of (3 — 4) for efa 1. Since the baryon den¬ 
sity has to be smaller than the number density of the AD field, e is at most unity even for 
large a# and fa. The amplitude of the flat direction decreases as time evolves due to the 
Hubble expansion and the B — L breaking effect is absent soon after the oscillation. Thus, 
the generated B — L asymmetry is conserved soon after the AD field starts to oscillate. 

Then, the oscillating AD field decays and dissipates into radiation [30] and the sphaleron 
effect converts the B — L asymmetry to the baryon asymmetry [12, 13]. Since the sphaleron 
process is in thermal equilibrium, the resultant baryon asymmetry is related to the B — L 
asymmetry as [31] 


nt = s”®- 1 - 


( 21 ) 


Assuming the absence of entropy production other than the reheating by inflaton decay, we 
can calculate the resulting baryon-to-entropy ratio Yj, as 


Y b = 




rib _ 8 n B -L 
s 23 s 

8.7 x HT n eg I 


^ 8 3 T RH n B _ L ^ 8 egT RH 

rh 23 4 p s osc 23 4 Hj 

Trh \ ( H: y 1 

4 x 10 3 GeV/ \4x 10 12 GeV J 



( 22 ) 


where ps (— 3i7|Mp L ) is the energy density of the inflaton S. This can be consistent with 
the observed baryon asymmetry of Y b obs ~ 8.7 x 10 -11 [1]. 

You can find differences from the conventional scenario of the Affleck-Dine baryogenesis. 
The Hubble parameter at the beginning of oscillation Hi is determined by the energy scale 
of inflation, not by the curvature of the potential for the flat direction (see Ref. [15], for 
example). This is because the flat direction starts to oscillate just after the end of inflation 
due to the positive Hubble-induced mass term, while in the conventional scenario it starts 
to oscillate at H(t) ~ m«^, where is the soft mass of the AD field. This allows us to 
consider a relatively large reheating temperature even if the initial VEV fa is as large as 
the Planck scale. In addition, the ellipticity parameter e, which describes the efficiency of 
baryogenesis, can be much smaller than unity when is smaller than the Planck scale. 
This is because the phase direction of the AD field is kicked by a higher-dimensional Kahler 
potential, which is highly suppressed for a small VEV of the AD field. However, as shown 
in the next subsection, the baryonic isocurvature constraint requires the initial VEV to be 
as large as the Planck scale. In that case, e is of order unity for an = 0(1). 

One might wonder if the energy density of the AD field dominates that of the Universe 
in the case that its initial VEV is as large as the Planck scale. This may be true in the case 












of conventional Affleck-Dine baryogenesis, in which the AD field starts to oscillate when 
the Hubble parameter decreases down to the soft mass of the flat direction. However, the 
energy density of the AD held never dominates the Universe in the above scenario because 
it decreases faster than that of radiation. Just after inflation, the AD held starts to oscillate 
around the origin due to the positive Hubble-induced mass term. Then, its number density 
decreases with time as a -3 due to the expansion of the Universe. This means that its energy 
density decreases as a~ 9 / 2 because its effective mass is of order the Hubble parameter, which 
decreases as a -3 / 2 . When the Hubble parameter decreases down to the mass of the AD held, 
that is, when H(t) ~ m^, its energy fraction to the total energy density is given as 


Pad 

r^j 

(mA 

Pad 

Ptot 

H—m^, 


Ptot 





Hr 


4 x 10 12 GeV 


(23) 


Thus, the energy density of AD held becomes negligible soon after inflation and the result 
of Eq. (22) is applicable to the case of ~ M Pl . 


3.3 Baryonic isocurvature perturbations 


Although the initial phase and radial values of the hat direction are almost constant over 
the whole range of the observable universe, they acquire quantum fluctuations like [9, 10, 
11, 14, 15] 


m\ ~ 


\S(pi 


Hj 

2tt \ 4 > i \ 

Hr 

_ 

~ 2vr‘ 


(24) 

(25) 


These fluctuations result in baryonic isocurvature perturbations because the produced baryon 
density is related to the initial phase 6i and VEV 0* (see Eqs. (19) and (20)). The baryonic 
isocurvature perturbation Sby is given by 


c SY b 
o bl = — ~ n 


Y, 


B 


^cot (nd) 56 + 



(26) 


Since we consider the case that the AD held starts to oscillate due to the positive Hubble- 
induced mass term, its radial direction also has quantum huctuations and contributes to the 
baryonic isocurvature perturbations. This leads to an additional factor in Eq. (26), which 
cannot be suppressed by the tuning of the initial phase 6 P Note that the VEV of the AD held 
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is smaller than the Planck scale due to the exponential factor in the supergravity potential. 
This leads to a lower bound on baryonic isocurvature perturbations like 


|Sb 7 | 




> 

r^j 


2.7 x 10" 7 x n 
2.7 x 10" 7 x n 


Hr 


4 x 10 12 GeV 
H r 

4 x 10 12 GeV 



? 


(27) 


where we assume 1 + cot (n6) ~ 1. 

Since the density perturbations of the cosmic microwave background are predominantly 
adiabatic [4, 5], the baryonic isocurvature perturbation is tightly constrained as [16, 15] 


\S hl \ <5.0 x 10" 5 . (28) 

Our scenario predicts the value below this constraint though it depends on the value of n 
and Hi (i.e., £). If one might consider a high-scale D-term inflation model, the resulting 
isocurvature perturbations can be as large as this lower bound and would be detected by 
future observations of CMB fluctuations. 


3.4 Comments on Q-ball formation 

In this subsection, we comment on Q-ball formation. If the potential of the AD field is 
shallower than the quadratic potential, its coherent oscillation is unstable and fragments 
into non-topological solitons, called Q-balls [32], The formation of Q-balls may change the 
scenario of Affleck-Dine baryogenesis [33, 34, 35, 36, 37, 38, 39]. For example, their decay can 
be another source of non-thermal production of DM [35, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 
50, 15, 51], or Q-balls can be a candidate for DM if they are stable [34, 52, 53, 54]. In the case 
considered in this paper, the AD field starts to oscillate by the positive Hubblc-induced mass 
term. When the beta function for the Hubble-induced mass of the AD field is positive, Q-ball 
does not form. The beta function has positive contributions from Yukawa interactions while 
it has negative ones from gauge interactions. The former positive contributions are roughly 
proportional to the squared masses of squarks and sleptons, and the latter negative ones 
are roughly proportional to the squared masses of gauginos. Here, since the Hubble-induced 
mass for gauginos is absent or 1-loop suppressed, the positive contributions from Yukawa 
interactions are usually dominant. Therefore, the beta function for the Hubblc-induced mass 
of the AD held is usually positive and Q-balls may not form in our scenario. However, if the 
AD held consists only of the hrst and second family squarks and/or sleptons, the positive 
contributions from Yukawa interactions are suppressed by small Yukawa couplings. In this 
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( 29 ) 


case, Q-balls might from. We estimate the typical charge 0 of Q-balls as 

Q ™ & {———' 

\ rn <j>,eff 


where m^^ff is the effective mass and f3 (~ 1CD 2 ) is a numerical factor obtained from 
simulations of Q-ball formation [36, 37, 55]. Here, we should substitute the Hubble-induced 
mass into the effective mass, and so the typical charge of Q-balls is at most 10 8 . Such small 
Q-balls soon evaporate into thermal plasma via interactions with the thermal plasma [58, 59] 
(see also Ref. [38]). 5 Therefore, the subsequent cosmological scenario and the calculation of 
the baryon asymmetry does not change. 

Even if Q-balls do not form just after the end of inflation, they may form at the time 
of H(t) ~ After that time, the potential of the AD field is dominated by its soft mass 
term. If the beta function of the soft mass is negative, the AD field becomes to fragment 
into Q-balls at that time. Since rib oc H(t)<p 2 (t ) oc a~ 3 oc H(t) 2 until the Hubble parameter 
decreases down to the soft mass, the amplitude of the AD field at H(t ) ~ m$ is given as 


ITj 


1/2 


^1 #(9— m 4> 

This implies that a typical charge of Q-balls is given as 


Q ~ /3 — ~ ^ 

™ 4 >. 




(30) 


(31) 


This is at most 10 18 for typical parameters. Such small Q-balls are evaporate into thermal 
plasma soon after they form. Even if Q-balls survive, they are so small as to decay into 
quarks before the BBN epoch. However, they usually decay after the electroweak phase 
transition [56, 57]. Since the sphaleron process is decoupled at that time, the AD field has 
to carry a nonzero baryon charge (not B — L) to generate the baryon asymmetry. In that 
case, the resulting baryon asymmetry is given by Eq. (22) without the factor of 8/23. 


4 Model for solution to the baryon-DM coincidence 
problem 

In this section, we propose a D-term inflation model which predicts an 0(1) ratio of baryon 
to DM density. We introduce a shift symmetry to ensure the flatness of the inflaton potential 

5 The evaporation is efficient during the inflaton oscillation era. In addition, since the energy per unit 

charge for these Q-balls is given by the Hubble parameter, their energy density decreases with time as a~ 9 ' 2 . 

Thus, the energy density of the Q-balls never dominate that of the Universe. 
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above the Planck scale. When there exists a small linear term in the Kahler potential, the 
inflaton decays mainly into gravitinos [18, 19, 20, 21, 22, 23, 24], The subsequent decay of 
those gravitinos is a source of non-thermal production of LSP DM and the resulting DM 
abundance is proportional to the reheating temperature and inversely proportional to the 
inflaton mass. Since the amount of the baryon asymmetry in Eq. (22) has similar parameter 
dependences, the baryon and DM densities are related with each other through the energy 
scale of inflation. 

In the next subsection, we propose the D-term inflation model. In Sec. 4.2, we investigate 
reheating processes of the D-term inflation model and then calculate the abundance of DM. 

4.1 Model 

Let us introduce a shift symmetry and an approximate Z 2 symmetry for the inflaton held 
S [17]. Under these symmetries, S transforms as S —* S + ict (cereal) and S —* —S, 
respectively. Then, the Kahler potential is written as 

K = c s (S + S*) + \ (S + S*f + |fM 2 + IVtI 2 , (32) 

where cs (<C 1) is an order parameter for the Z 2 symmetry breaking effect. The superpo¬ 
tential of Eq. (2) explicitly breaks the shift symmetry, which is required to ensure a graceful 
exit. Otherwise the inflaton stays at a certain VEV because it has a exactly hat potential. 
In this model, we should replace \S\ 2 with (S + S*) 2 /2 for the calculations in Sec. 3.1, though 
the results are unchanged. 

There is an advantage to impose the shift symmetry to the inflaton. In order to obtain a 
sufficiently large e-folding number, say, TV* ^ 60, the initial VEV of the inflaton S has to be 
as large as A r *M P] ~ 0.5Mpi. which is of order the Planck scale. This implies that the 
Planck-scale physics may affect the potential of the inflaton and spoil its flatness. However, 
the shift symmetry ensures the flatness of the inflaton potential above the Planck scale. 

If Z 2 symmetry is exact, some MSSM particles have to carry odd Z 2 charge and interact 
with the field S [17] for the inflaton to decay. In this case the LSP DM abundance is given 
by the usual thermal relic density. Here, we introduce Z 2 breaking terms in the Kahler 
potential so that the field S efficiently decays into gravitinos [23, 24], whose decay is a 
source of non-thermal production of LSP DM. 

We assume that the mass of gravitinos is of order 10 2 ~ 3 TeV so that gravitinos decay 
into radiation before the BBN epoch. Otherwise the decay of gravitinos spoils the success 
of the BBN, or their energy density overdoses the Universe if they are stable. Such a heavy 
gravitino is well motivated in a class of SUSY models with a split spectrum [60, 61, 62, 63, 
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64, 65, 66]. In these models, the masses of gravitino as well as squarks and sleptons are of 
order (or larger than) 10 2-3 TeV while those of gauginos are of order 1 TeV. This hierarchy 
can be realized when gauginos acquire one-loop suppressed soft masses through the anomaly 
mediated SUSY breaking effect [67, 68]. In that case, the mass of wino and gravitino is 
related with each other such as 

ma, = (m 3 / 2 + L) ~ 3 x 1CT 3 (m 3/2 + L) . (33) 


The factor L is the Higgsino threshold corrections and is calculated as [68, 69] 


L = ijlh sin 2/3- 


m. 


log 


|/i 


Hi 


m. 


(34) 


M' - m A 

where mA is the mass of the heavy Higgs bosons, /jlh is the SUSY mass of the higgsinos, 
and tan /3 is the ratio of the VEV of H u and H^. When /i# is of order the gravitino mass, 
the Higgsino threshold correction is important and the wino mass is ~ 10 _3 m3/2. Note that 
neutral higgsino can also be the LSP when /i# is sufficiently small. The following discussion 
does not rely on the detailed properties of the LSP except for its mass. Hereafter, we assume 
that the mass of gravitino is O(10 2 ~ 3 ) TeV and that of the LSP is O(10 2 ~ 3 ) GeV. 

In order to calculate the gravitino production rate from the inflaton decay, we need to 
specify the SUSY breaking sector. We introduce a Polonyi held z, which breaks SUSY in 
low energy scale, and consider a simple extension of the Polonyi model given as 

|4 


K = |H 2 - 


\z\ 
A2~’ 


W = fz + W„, 


(35) 

(36) 


where A is a cutoff scale, /i is the SUSY breaking scale, and Wo is a constant term which makes 
the cosmological constant (almost) zero in the present Universe. This can be achieved by the 
O’Raifeartaigh model after integrating out relatively heavy particles [70] or by dynamical 
SUSY breaking models, including the IYIT model [71, 72]. The important parameters are 
calculated as 


/l 2 

- V3m 3 / 2 M Ph 

(37) 

m 2 z 

12 m 2 /2 

- 

(38) 

z )o 

~ 2^3 ( m3/2 ) M P1 , 

\ m z J 

(39) 


where m z is the mass of z and ( z ) is its VEV at the low energy vacuum. Since the Hubble- 
induced mass is absent during D-term inflation, massless scalar fields cannot be stabilized 
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at the origin. This implies that if the mass of the Polonyi field is much smaller than the 
Hubble parameter, it obtains a VEV as large as the Planck scale during inflation. In this 
section, we consider the case that the mass of the Polonyi held is larger than Hj and it stays 
at the origin of the potential during inflation, while we consider the case of relatively light 
Polonyi in Appendix. The conditions of m z > Hj and A 0 // in the effective theory leads to 
the lower bound on the gravitino mass [73]: 


^3/2 ~ 




V3 H] 

12Mpi 

10 3 TeV 


Hr 


4 x 10 12 GeV 


(40) 


Thus, the heavy gravitino is favoured in D-terrn inflation to stabilize the VEV of the Polonyi 
held. 


4.2 Reheating process 

In this subsection, we investigate the reheating process and calculate the DM abundance in 
the model introduced in the previous subsection. 

As explained in Sec. 2, the held ip + decays into the MSSM helds much faster than the 
inhaton S, so that the reheating temperature of the Universe is determined by the relatively 
late-time decay of the inhaton S [8]. After the held 0 + decays completely, the effective 
superpotential can be rewritten as 

IU (inf) = m s S^_, (41) 

m s = (42) 

where ms is the effective mass of the helds S and 0_. This superpotential is equivalent to 
the one in the model of chaotic inflation proposed in Ref. [17], except for the value of ms- 
The supergravity effects induce a soft SUSY breaking B-terrn of &m 3 / 2 msS'0_, where b 
is an 0(1) constant. This implies that they maximally mix with each other and form mass 
eigenstates 

$± = V (s ± </i) , (43) 

around the potential minimum [19, 23]. Therefore, when the time scale of inhaton decay 
T" 1 is longer than that of the mixing effect m^, we have to consider the decay of $-t- 
to investigate the reheating process. Since we consider a heavy gravitino and a reheating 
temperature of order 10 3 GeV (see Eq. (22)), the mixing effect is indeed relevant. 
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The Z 2 breaking term in the Kahler potential results in the decay of the field $_)_ through 
supergravity effects [24], 6 First, let us focus on the top Yukawa interaction in the MSSM 
sector: 

kb (top) = y t Q 3 H u u c 3 , (44) 

where yt is the top Yukawa coupling constant, and Q 3 , H u , and u 3 are the chiral supermulti- 
plets of the MSSM sector. The relevant interaction terms between 0_ and the MSSM fields 
are given by 

U = i ^-K s W ( - top) W* s + c.c. + ..., 

IVlpi 

= ytn ^2 <S {Q:i H v.u-i) + C.c. + ..., (45) 

m pi 

where the dots represents the other irrelevant terms. Since the fields $- 1 - consist of 

0_ as Eq. (43), they decay into the MSSM scalar fields through this interaction. They also 
decay into the MSSM fermion fields, which equally contributes to the $- 1 - decay [24], Thus, 
the partial decay rate of $- 1 - into the MSSM fields is given as 

Q 7T7 ^ 

Tmssm (»± -> MSSM) = ^ |»| 2 (46) 

where we use Ks = cs■ Since we consider the gaugino mass (rrij) much smaller than the 
gravitino mass, the decay rates of S into gauge fields are suppressed by a factor of {rrig/m 3 / 2 ) 2 
and can be neglected [22], 

Next, let us consider the decay of <f>± into gravitinos [18, 19, 20, 21, 22, 23, 24], We 
follow the discussion presented in Ref. [74], When the field S has a nonzero VEV, the field 
0_ mixes with the SUSY breaking field z and can decays into goldstino, (i.e., longitudinal 
component of gravitino). This is because the supergravity effects induce mixing terms such 
as 


V = W s (K S W)* + KglW s W* + c.c. + ... 

= m s dF z iJ)_z* + c.c. + ..., (47) 

d = (K s )-(K S z2 ), (48) 

where the dots ”...” represent the other irrelevant terms. The second term in d is relevant 
when there is a term like (S + S*) \z\ 2 in the Kahler potential, whose coefficient is of order 

6 Since the field 0_ has a small VEV of order cm 3 / 2 /ms due to the supergravity effects, it can decay into 
the MSSM fields through the D-term potential. However, we confirm that its partial decay rate is irrelevant 
due to the suppression factor coming from its small VEV and can be neglected. 
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eg ■ Thus, the fields i\)_ and z mix with each other and the mixing angle is given by 


6 ~ 


d 


F z m s 


m 


2 ’ 
z 


( 49 ) 


where we use m z ms- Since the fields <h-t consist of , 0_, they mix with z and the mixing 
angle is given by 9/-\[2. Since the SUSY breaking field z has an operator of 


£ = —2 ^zz^zt + h.c., (50) 

A 2 

it decays into goldstino z. Together with the mixings between $_)_ and z, the fields decays 
into goldstino through this operator. The partial decay rate of the field into goldstino is 
therefore calculated as [75] 

r 5 ($±-^i) ~ 


From Eqs. (46) and (51), the total decay rate Tg and the branching ratio of the decay 
of $- 1 - into gravitinos S 3/2 are given by 


1 

32tt 


d 2 m' 


_e_ 

7! 

3 

s 


ml 


J m Si 


64v tM^' 


(51) 


Tg 

Br 3/2 


Tmssm ($± MSSM) + T s ($± ->■ zz ), 
d 2 

d 2 + 3 \y t \ 2 c 2 / (in 2 ) 


(52) 

(53) 


Since d/c = 0(1) and yt = 0(1), the branching ratio is almost unity. This means that the 
energy density of the Universe is dominated by that of the gravitinos after the fields 4>-t (i.e., 
the inflaton S) decay completely. ' 

Since the fields $- 1 - are much heavier than gravitino, the produced gravitinos are highly 
relativistic. The Lorentz factor for the gravitinos at a time i/~ 1 (f) is given as 


l{t) 


m s \ 

m 3 / 2 ) 


H(t) 

Tg 


1/2 


+ 1 




ms 

m 3/2 



(54) 


The gravitinos decay into MSSM particles with a rate of 

, 7 . ,r 5 

- 7 _1 W 


1 4 


3/2 


487T 771^2 Mpj ’ 


(55) 


7 Note that since we consider relatively low reheating temperature ~ 10 3 4 GeV, we can neglect the 
thermal production of gravitinos [76, 77]. 


16 














where the summation is taken for all MSSM particles X,. Since we consider a SUSY model 
with relatively light gauginos and relatively heavy squark and sleptons, we can roughly 
estimate the numerator as 24m3^ 2 . This implies that the gravitino decays into radiation at 
the temperature 



where g * (~ 10.75) is the effective number of degrees of freedom at the decay time. We 
require that the mass of gravitino is of order 10 2 “ 3 TeV or larger so that its decay completes 
before the BBN epoch, that is, T 3/2 > 1 MeV. Otherwise the decay particles interact with 
the light elements and spoil the success of the BBN [78, 79, 80, 81, 82], The gravitino decay 
temperature T 3 / 2 is much smaller than the mass of the LSP, so that the decay of gravitino is 
a source of its nonthermal production. Since the energy density of the Universe is dominated 
by that of gravitino before they decay, the thermal relic density of the LSP is diluted by the 
entropy production from the gravitino decay. Therefore, the LSP abundance is determined 
by the nonthermal production from the gravitino decay. The produced number density of 
the LSPs is equal to that of the gravitinos due to the R-parity conservation. Note that the 
annihilation of the produced LSP is usually inefficient in such a low temperature. 

The Lorentz factor of the gravitino is of order 10 3 for the reference parameters shown in 
Eq. (56). This implies that the scale factor of the Universe continues to decrease as a -1 from 
the time of reheating by the decay of $- 1 -. Although the LSPs are relativistic at the time 
they are produced from gravitino decay, they lose their energy through interactions with the 
thermal plasma and soon become to non-relativistic particles [83, 84, 85]. Therefore, the 
LSP DM is cold even though they are produced non-thermally in this scenario. 


4.3 DM density and baryon-DM coincidence 

Let us summarize the scenario of non-thermal production of DM. First, the inflaton S (or 
$±) decays into gravitinos as well as the MSSM particles at H(t) — Ts- Then the energy 
density of the Universe is dominated by the relativistic gravitinos and decreases as a -4 . The 
gravitinos decay into the MSSM particles just before the epoch of the BBN and the LSP 
DM is produced non-thermally. Since the thermal relic density of the LSP is diluted by the 
entropy production of gravitino decay, its abundance is determined by the gravitino decay. 
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Thus, we can estimate the resultant DM abundance as 

™LSP 


Edm = 


s 

n 3/2 


h =r 


3/2 


3T 3 /2 Tl 3 /2 


4 

3 ^ 3/2 


P3/2 

r s 


(ft) 


h= r 3/2 
1/2 


^3/2 


P3/2 


n=r s 


3T^ e H ff) 2Br 3/2 n s 


Ps 


H=r s 


3T, 


(eff) 


RH 


2m,c 


(57) 


where we have used Br 3 /2 — 1 in the last line. We define the effective reheating temperature 
by Eq. (10) with the replacement of g*(T RH ) —> g*(T 3 / 2 ) as 

1/4 


T, 


(eff) 


RH 


90 


,g*(T 3/2 )n 2 
~ 1.5 x 10 3 GeV 


\/IWpi 


m s 


\ 3/2 




d 


10 


-10 


(58) 


5 x 10 15 GeV^ 

The reheating temperature is adjusted by the Z 2 symmetry order parameter d to obtain a 
desirable amount of baryon asymmetry from Eq. (22) or DM from Eq. (57). 

Here we take into account the baryon asymmetry generated by the Affieck-Dine baryo- 
genesis. Once we replace the reheating temperature T RH with the effective one T^p defined 
by Eq. (58), the resulting baryon asymmetry is still given by Eq. (22) even in this scenario. 8 
Combining Eqs. (22) and (57), we obtain the following simple relation for the baryon-to-DM 
ratio: 

Qb 4 m r , 




DM 


" p TYlS 

-TT’ 

69 m L sp H] 


(59) 


where we assume 0* ~ Mp\. 
m p ~ 0.938 GeV, we obtain 


Substituting benchmark parameters and the proton mass 


Qh 


Q 


DM 


~ 0.22eg 


"Ilsp \ 


-l 


m s 


400 GeV/ I 6.6 x 10 15 GeV \ 4 x 10 12 GeV 


Hr 


-l 


8 Note that the inflaton also decays into the MSSM fields with a branching ~ 10 - d -2 \ so that there 
exists significant thermal plasma after the decay of inflaton. Thus, the sphaleron effect proceeds fast enough 
to convert the B — L asymmetry to the baryon asymmetry even if the dominant component of the Universe 
is gravitino at that time. 
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(60) 


_ O.UeqXg” 1 (- 

V 




-i 




m LS p 

400 GeV J V 6.6 x 10 15 GeV 


-l 


which is naturally of order unity and is consistent with the observed value of ~ 

0.2 [1]. The scenario naturally explains the coincidence of their energy density, known as the 
baryon-DM coincidence problem. This is because both of them are related to the energy scale 
of inflation. The amount of baryon asymmetry is proportional to the reheating temperature 
of the Universe and inversely proportional to the Hubble parameter during inflation. That 
of DM is proportional to the reheating temperature and inversely proportional to the mass 
of inflaton. Since the Hubble parameter (77/ ~ g^/Mp{) and the inflaton mass (ms = A-\/£) 
is related to each other, the resulting baryon and DM density is naturally of order unity. 
Interestingly, the electroweak scale DM mass (O(10 2 ) GeV) comes from the fact that the 
GUT scale, which is expected to be, is two orders of magnitude less than the Planck 
scale. 

Although the result has an 0(1) uncertainty coming mainly from A and £, the LSP with 
mass of O(10 2 ~ 3 ) GeV is favoured in our scenario. If the LSP DM is mostly wino or higgsino, 
the indirect detection experiments of DM puts lower bounds on DM mass. The wino DM 
with mo, < 390 GeV and 2.14 TeV < m& < 2.53 TeV is excluded [86], while the higgsino 
DM with m^ < 160 GeV is excluded [87]. The future indirect detection experiments can 
detect the wino DM with m& < 1.0 TeV and 1.66 TeV < < 2.77 TeV [86]. 


5 Summary 

We have considered the Affleck-Dine baryogenesis assuming a positive Hubble-induced mass 
term after D-terrn inflation. We have calculated the baryon asymmetry and found that 
the result is consistent with the observed abundance of baryon asymmetry for reheating 
temperature Trh ~ 10 3 GeV. There are some differences from the conventional scenarios of 
the Affleck-Dine baryogenesis where the sign of the Hubble-induced mass term is negative. 
First, since the AD field starts to oscillate just after the end of inflation, the resulting baryon 
asymmetry is inversely proportional to the inflation scale 77/. This allows us to consider a 
relatively large reheating temperature even if the initial VEV <f>i is as large as the Planck 
scale. In addition, the ellipticity parameter e, which describes the efficiency of baryogenesis, 
can be much smaller than unity when fa is smaller than the Planck scale. This is because 
the phase direction of the AD field is kicked by a higher-dimensional Kahler potential, which 
effect is highly suppressed for a small VEV of the AD field. Since the radial direction as 
well as the phase one has quantum fluctuations during D-terrn inflation, the Affleck-Dine 
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baryogenesis predicts nonzero baryonic isocurvature perturbations. They would be detected 
by future CMB observations if one considers a high-scale D-term inflation model. 

We also proposed a D-term inflation model with a shift symmetry in the imaginary 
direction of the inflaton superfield and a small linear term in the Kaliler potential. We 
consider the case that the mass of gravitino is O(10 2 ” 3 ) TeV and the mass of the LSP 
is O(10 2 - 3 ) GeV, which is naturally realized in anomaly mediated SUSY breaking mod¬ 
els. In this model, the inflaton decays mainly into gravitinos through the supergravity 
effects [18, 19, 20, 21, 22, 23, 24] and the subsequent decay of gravitinos is a source of non- 
thermal production of DM. Together with estimation of baryon asymmetry generated from 
the Afileck-Dine mechanism, the resulting DM density gives an 0(1) baryon-to-DM ratio. 
This is because both of them are related to the energy scale of inflation. The amount of 
baryon asymmetry is proportional to the reheating temperature of the Universe and inversely 
proportional to the Hubble parameter during inflation. That of DM is proportional to the 
reheating temperature and inversely proportional to the mass of inflaton. Since the Hubble 
parameter and the mass is related to each other, the resulting baryon and DM density is 
naturally of order unity. We predict that the LSP mass is two orders of magnitude larger 
than the proton mass, which comes from the fact that the GUT scale is two orders of mag¬ 
nitude less than the Planck scale. When the LSP is mostly wino or higgsino, it would be 
detected by future indirect detection experiments of DM. 
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A Another solution for baryon-DM coincidence prob¬ 
lem 


In this appendix, we consider the case that the mass of the Polonyi field is less than the 
Hubble parameter Hi and it obtains a nonzero VEV during inflation. Even if the origin of its 
potential is a symmetry-enhanced point, the Polonyi field cannot be stabilized at the origin 
because of the absence of the Hubble-induced mass during D-term inflation. Although one 
might wonder if this is an obstacle known as the Polonyi problem, we show that its decay 
can explain the amount of DM once we allow an 10% fine-tuning for the initial VEV of the 
Polonyi field. 

During D-term inflation, the Polonyi field obtains a VEV denoted as z t , which might be 
as large as the Planck scale. Just after inflation ends, the Polonyi field starts to oscillate 
around the origin due to the positive Hubble-induced mass term as the AD field considered 
in Sec. 3. Then, its number density decreases with time as a~ 3 due to the expansion of 
the Universe. This means that its energy density decreases as a -9 / 2 because its Hubble- 
induced mass is of order the Hubble parameter, which decreases as a~ 3//2 . After the Hubble 
parameter decreases down to its low-energy mass of the Polonyi field, that is, after the time 
of H(t) ~ m z , the Hubble-induced mass term can be neglected. Then its mass and VEV are 
given by Eqs (38) and (39), respectively. Note that the minimum of the potential is usually 
much smaller than the amplitude of the Polonyi field at that time because its amplitude 
decreases only as a~ 3 / 4 until H(t) ~ m z : 


AJ \H(t)=m z 



1/2 

Zj » (z)„. 


(61) 


This means that its number density is not affected and continues to decreases with time as 
a" 3 . 

The Polonyi field decays mainly into gravitinos with a rate of 


r Z {z —$■ 2 - 03 / 2 ) — 


1 ml 

967T 


(62) 


Since we consider relatively low reheating temperature to realize the Affleck-Dine baryoge- 
nesis, we can neglect the thermal production of gravitinos. Thus, the gravitino abundance 
is determined by the number density of the Polonyi held. We require that the mass of grav¬ 
itino is of order 10 2 ~ 3 TeV so that it decays into the MSSM particles before the BBN epoch. 
Otherwise the decay particles interact with the light elements and spoil the success of the 
BBN [78, 79, 80, 81, 82], The decay of gravitino is a source of non-thermal production of 
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LSP DM. We assume that the thermal relic of the LSP is much smaller than the observed 
DM abundance. This can be achieved for the case of wino-like LSP with a mass much less 
than 3 TeV [88, 89] or higgsino-like LSP with a mass much less than 1 TeV. The wino-like 
LSP is well motivated in models of anomaly mediated SUSY breaking, where the gravitino 
mass is naturally as large as 0(100) TeV as we required. 

In summary, the LSP DM is nou-thermally produced from the decay of gravitino, which is 
generated from the decay of the Polonyi field. The Polonyi held is generated by its coherent 
oscillation just after the end of inflation. We thus obtain the following DM abundance: 


Ydm = 












riLSP 

s 

^ 3/2 


S 


ip 3/2 decay 


3Trh%/2 


4 pi 


3Trh7U 


2 Pi 

h 

3Trh^z 


2 pi 

H 

TrhU 2 



2H I M^ 


(63) 


where we use n z ~ Hjzf at H = Hi in the last line. 

Together with estimation of baryon asymmetry generated by the Affleck-Dine mechanism 
calculated in Sec. 3.2, the above result implies that the baryon-to-DM ratio is given by the 
following simple relation: 


Qb ^ 4 ^ mp f <k_ 
Ddm 23 mLSP \ z i 


(64) 


Thus, we can explain the observed baryon-to-DM ratio when ~ 0.1. One might expect 
that the natural values of their initial VEVs are of order the Planck scale. In that case, the 
result requires a 10% fine-tuning for the value of Z{. 
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